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DRAG FORCES IN VELOCITY GRADIENT FLOW* 


7 Frank D. Masch, A. M. ASCE! and Walter L. Moore, M. ASCE2 


SYNOPSIS 


n exploratory investigation was made of the drag coefficient for a circular 
ider as influenced by a velocity gradient along its axis. Someof the three- 
nsional aspects of the flow mechanics were described qualitatively. Ex- 
ments indicated that the local drag coefficient varied significantly along the 
hof the cylinder. For velocity gradients inthe range of engineering signi- 
ice, local drag coefficients may be affected by as much as 40%. 


INTRODUCTION 


he evaluation of the hydrodynamic forces on a solid body in a flowing fluid 
primary importance in many engineering problems, In afew simple cases, 
> forces can be evaluated by analytical methods alone. In most practical 
lems, it isnecessary to resort to a combination of analysis and experiment. 
<perimentally determined drag coefficient is commonly used in estimating 
forces. 

<perimental determinations of drag coefficients for various shapes usually 
been made in a stream with a nearly constant velocity. However, many 
lems of engineering importance arise in which the velocity of the fluid 
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stream is not constant. For example, velocity variations will be present if 
fluid movement is the result of wave motion, or if the fluid motion is retar 
by the resistance along a boundary. More specific problems, in which a vel 
ity gradient may be significant, are; wave forces’on piers or other structur 
wind forces on buildings or stacks, resistance due to piers in streamsa 
channels, forces on sediment particles, and dragforces on surfaces irreguli 
ities in a boundary layer, ® 

Comparatively little effort has been devoted to the study of drag forces i 
velocity gradient. S. F, Hoerner, (1)3 has treated two aspects of the proble 
He summarizes the work of Tillman, Wiegard, and others, on the drag of p! 
tuberances ina turbulent boundary layer. For these surface irregulariti 
whose height is less than the boundary layer thickness and whose widthis la} 
in relation to their height, the indications are that a satisfactory estimate 
the drag force can be made assuming the drag coefficient tobe the same as - 
flow with no velocity gradient. An effective drag coefficient based on the vy 
locity at the top of the protuberance can be computed by simply summing t 
forces at successive layers accounting only for the variation of the velocit 
That is, at each layer, the drag coefficient is considered to be the same as f 
two-dimensional flow. Hoerner also presents approximate data on the veloc: 
profiles near the ground for wind and on drag coefficients for several buildi 
shapes. 

A number of papers are available, dealing with the motion of solid bodi 
in rotating fluids. H. B, Squire’s (2) has dealt with several aspects of tl 
problem and refers to a group of papers on this subject presented in the peri 
from 1915-25, | 

W. R. Hawthorne (6) made a theoretical and experimental study of the fl 
about struts and airfoils in aflow with a spanwise velocity gradient. His ani 
ysis of vortex development is of special interest, but as he points out it is] 
applicable to blunt-nosed bodies. He assumes that the induced spanwise veloc 
components are small and he does not account for any effect in the wake, F 
these reasons his analysis, although basedon similar concepts of the mechal 
of the flow, is not applicable to the present investigation. 


FORMULATION OF PROBLEM 7 


The variables normally considered in drag force studies canbe grouped 
a dimensionless equation of the form 


FD 


Cp = (R, form) 


Ap V2/2 


in which Cp is the drag coefficient, Fp the drag force, A thecross-secti 
area of the object taken perpendicular to the flow, p the mass density of 
fluid, V the velocity of the approaching fluid, and R the Reynolds Number, 
This equation represents the average drag coefficient as a function of 
Reynolds Number and the shape or form of the object under consideration, 
application of this equation to two-dimensional flow problems normally assu 
a uniform velocity distribution approaching the object in question, 


3 Numerals in parentheses, thus (1), 
graphy. 


refer to corresponding items in the Bi 
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his investigation involved the introduction of additional variables in Eq. 1 
rmit an evaluation of the drag forces ona circular cylinder with its axis 
al to the approaching flow as affected by a linear velocity gradient in the 
ction of its axis, 

dimensional analysis enables a velocity gradient parameter to be formu- 
| into the drag coefficient equation as 


Fp 


he geometrical aspects are represented by the form and the ratios of cyl- 
‘diameter tocylinder length, D/Yo. If theform is restricted to a circular 
der and the ratio, D/ Yo, is held essentially constant, these two terms will 
ffect the drag coefficient. AV D 

he velocity gradient parameter, Wv” is a dimensionless combination of 


bles inwhich AV /AY isthe slope of the velocity profile and is is the ratio 
e diameter of the cylinder to the velocity measured at the midpoint of the 
der. If the velocity distributionis uniform, the slope, AV/AY is zero, thus 
ing a to vanish leaving the drag coefficient a function of the Reynolds 
ber only. 

1e Reynolds Number may be formed from the velocity at the mid-point, as 
mined from the velocity profile measurements, the diameter of the cyl- 
', and thekinematic viscosity of thefluid. For a rangeof Reynolds numbers 
2en 10,000 and 100,000, there is practically no variation of the average drag 
icient with the Reynolds Number. In this range of Reynolds numbers, the 
“forces on the cylinder are very small in relationto the pressure forces, 
he total drag can be considered as being primarily form drag. 

ider the limitations previously described, the form, D/Yo, and R may be 
nated from the functional relationship of Eq. 2, leaving the drag coeffi- 
a function of the velocity gradient parameter only. 

3 an aid in analyzing the mechanics of this problem, a local drag coeffi- 
was formulated similar to the average drag coefficient, that is, 


ich f is the local drag force, AY D is the projected area of an element of 
linder, and v is the local velocity of the fluid approaching the element. 
cal force, f , may be determined from an integration of the pressure dis- 
ion around the cylindrical element, or 


27 
f = ile PlPC68 0d Oe at ern She (4) 


ich p is the pressure intensity at any pt., r the cylinder radius, and 6 the 
ar position on the cylinder measured from the upstream element of the 
ler. 

mbination of this local force with the corresponding measured local ve- 
will give the local drag coefficient. This computation, repeated at several 
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elements along the longitudinal axis of the cylinder, will give the variatio 
the local drag coefficient along the cylinder. 

A consideration of the mechanics of flow around the cylinder will indi 
that the velocity gradient will induce flow components along the axis of the « 
inder. Along the upstream element of the cylinder (the stagnation line) , 
stagnation pressure will be greater at the endof the cylinder, where the vel 
ity is high, than at the end where the velocity is low. This will produce a pr 
sure gradient along the axis of the cylinder andinduce a flow along the stag 
tion line toward the low-velocity end of the cylinder. In a similar way, the 
duced pressure in the wake will be affected by the local approach velocity. 
the high-velocity endof the cylinder, the pressure in the wake will beless 1 
at the low-velocity end, thus, inducing a longitudinal flow toward the high-vel 
ity end of the cylinder, in the wake zone. 

The effect of these longitudinal flows on the local drag coefficient at vari 
locations along the cylinder may be anticipated asfollows. At the high-velo 
end of the cylinder the longitudinal flow along the stagnation line would be 
verting fluid away and, thus, reducing the pressure intensity on the upstrt 
side of the cylinder. At the same time, the longitudinal flow in the wake z 
would be forcing fluid toward the high-velocity endand thereby raising the pr 
sure intensities on the downstream side of the cylinder. Both of these eff 
wouldact to reduce the drag force onthe cylinder and, hence, to reduce the kl 
drag coefficient at the high-velocity endof the cylinder. It might be anticip: 
that the opposite effect would be true at the low-velocity end of the cylin 
and the local drag coefficient would be high. However, it is not necessa 
true that this value will exceed the one for normal two-dimensional flow. 

This rather crude description of the secondary flows may help to exp 
some of the essential features of theflow mechanics, but itappears that a m 
adequate description can be developedinterms of vortex generation. Some ¢ 
cepts bearing on this approach are described by Hawthorne (6) and by Sq 
and Winter (10) but more study is needed to developthem for application to 
problem of this investigation. 


EXPERIMENTAL INVESTIGATION 


Experiments were made to explore the variation of the local drag coeffic 
along a cylinder, A 1-in. diam hard copper tube was mounted vertically at 
center of a 12-in. wide glass-walled flume in which a flow of water was m: 
tained at adepth of 0.8 ft. The flow rate was adjusted to give Reynolds Numk 
at mid-depthin the range between 18,500 and 22,500,a region where the nor 
two-dimensional drag coefficientis constant. With the Reynolds Number in 
range, and the form and D/Yo held constant, these factors were eliminate 
a source of variation in the measured drag coefficients. Fig. 1 is a schem 
diagram of the test arrangement. 7 

Since the flow had a free surface it might appear that the drag coeffici 
would be influenced by the free surface and hence the Froude Number sl 
be considered. This would of course be true if the total drag coefficient 
under consideration because the drag produced by the wave at the surface 
be a very important part of the total drag. The investigation reported 
deals with the local drag coefficient as determined from pressure measure 
on the part of the cylinder that was completely submerged. It is believe 
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free surface had a negligible effect onthe pressures measuredat the cylin- 
calsurface except for points near thefree surface. It is, therefore, believed 
t the free surface had no effect on the results reported here. 

Apparatus .— Control of the velocity distribution was obtainedwith a variable 
tter gate. Variation of the gate opening from top to bottom permitted the 
istance to be varied from a maximum near the bottom to a minimum at the 
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FIG, 1.—DIAGRAM OF TEST ARRANGEMENT 


ace. The effectiveness of the gate in producing a linear velocity variation 
enhanced by the inclusion of a serrated screen. 

the cylinder was mounted to an aluminum plate suspended by four thin alu- 
um straps from arigid frame. (Fig. 2) The lower plate was free to deflect 
orce was applied to the cylinder. Deflection of the lower plate introduced 
ural stresses in the aluminum straps, the strain from which was picked up 
system of electrical strain gages, thereby enabling total force measure- 
ts to be obtained. 
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whe cylinder was attached to the lower plate by a sleeve fitted into a collar 
d to the lower plate. It was possible to rotate the cylinder through 360° 
it its longitudinal axis. Located along a vertical line on the cylinder were 
mezometers, each of which was connected internally toa 1/8 in. copper 
. The tubes were brought out through the top of the cylinder andconnected 
manometer board from which the local pressure was measured for each 
ometer. Fig. 2 shows the cylinder and the amounting assembly. 
xperimental procedure.—The experimental procedure involved measure- 
ts of the velocity distribution, local pressure distribution, and total drag 
e. The longitudinal velocity profiles were measured with a Prandtl pitot 
Several transverse profiles were also obtained to insure that side wall 
sts were negligible. Pressure distributions around the cylinder were ob- 
sd with the piezometers by rotating the cylinder through 360° in increments 
)”, Total drag force measurements were taken several times during each 
All of these measurements were taken at four different values of the ve- 
y gradient parameter. Total drag measurements were made as a check on 
ntegrated pressure measurements. The measured total drag included, of 
se, the wave drag at the free surface. An approximate evaluation of the 
> drag was made, based on the profile of the wave surface around the cyl- 
r. The sum of the estimated wave drag and the pressure drag obtained 
| integration of the pressure measurements agreed satisfactorily with the 
drag force as measured by the electrical strain gages. 
he evaluation of the local drag coefficient at various points along the cylinder 
jetermined from the pressure distributions and the measured velocity pro- 
- To illustrate the method of determining the local drag coefficient, consider 
ction of the cylinder. The pressure distribution as determined from the 
ymeter at the section was integrated tofind theforce at the particular sec- 
As the functional relation between the pressure and the angle couldnot be 
essed mathematically, Eq. 4was evaluated by graphical integration. Typical 
e curves used for this integration are those of Fig. 3. This plot includes a 
e for each piezometer. 
aving completed the integration of the pressure distribution on a section, 
orresponding local velocity was taken from the plot of the velocity profile, 
4, The local drag coefficients were computed at various positions along 
ylinder for different values of the velocity gradient parameter inthe man- 
lescribed above. 


RESULTS OF INVESTIGATION 


eal drag coefficient.— The velocity gradient very definitely affects the local 
coefficient, The variation that exists may best be seen in Fig. 5. Here the 
drag coefficient is plotted against the distance from bottom in feet for 
‘ent values of the gradient parameter. If the velocity distributions were 
rm, there would be no variation in the local drag coefficient as the velocity 
ressure distribution would be constant along the cylinder. This curve in- 
28 that the local drag coefficient does vary with position and also with the 
of the slope of the velocity profile. 

msider first the variation with position along the cylinder. Fig. 5 indicates 
he local drag coefficient is a maximum at the bottom of the cylinder and a 
num near thetop. With reference tothe velocity gradient, this corresponds 
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_high coefficient at the low-velocity end and a low coefficient at the high- 
xcity end. This variation inthe drag coefficient was,as anticipated and pre- 
usly described, in terms of the secondary flows. 

There is, also, a variation of the local drag coefficient with the steepness 
he velocity gradient. From the curves in Fig. 5, it can be seen that as the 
be parameter increases from 0.058 to0,117, the variation of the local drag co- 
cient becomes more pronounced. This follows from the fact that the greater 
velocity gradient, the greater is the pressure variation in the wake zone, 
the greater is the upward flow in the wake zone. 

To illustrate the upward flow that exists in the wake, a dye was injected 
ugh one of the copper tubes andallowed to issuefrom a piezometer located 
he wake zone. The result of this qualitative demonstration is presented in 
. 6. The dye jet may be seen shooting out from the piezometer for a dis- 
se of about 1/4 of the cylinder diameter and then turning upward where it 
s the back of the cylinder. Although there is some subsequent diffusion of 


FIG. 6.—UPWARD FLOW IN WAKE. 


dye into the surrounding fluid as it is swept downstream, the existence of 
upward flow is emphasized by the fact that only a trace of the dye appears 
Ww the piezometer opening for the downstream region. 

verage drag coefficient.—In general, it may be concluded that the average 
- coefficient for a cylinder witha longitudinal velocity gradient is less than 
sonventional value of 1.2 for the chosen range of Reynolds Numbers. The 
lary reason for this is that the flow isnot two-dimensional. The presence of 
ongitudinal flow in the wake zone very definitely indicates three-dimensional 

characteristics. 
ne values of the average drag coefficient agree quite well with approx- 


e values given for cylinders of finite length. For values of 5< = e< 20, the 


we drag coefficient varies from 0.74 to 0.90. (3) The results R, this in- 
igation more nearly approach these values at the steep velocity gradients, 
ver, the causes are not directly related. 


SIGNIFICANCE FOR ENGINEERING PROBLEMS 


finvestigation has shown that the longitudinal velocity gradient can effect 
ue of the drag coefficient and that a longitudinal flow is induced in the 
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wake of the cylinder which appears to be related to the effect. It is of inte 
to compare the values of the velocity gradient parameter used in the investiga 
with values that might be expected in some engineering problems. In the 
perimental investigation the parameter AV/AY D/V had the following vah 
0.117, 0.077, 0.061, and 0.058. 

For comparison, some values of the velocity gradient parameter were ¢ 
culated for waves based on laboratory measurements of particle velocitie 
waves by T. A. Marlowe. (4) Ifitis assumed thata pier diameter is 0.05 til 
the water depth, the wave length is 2.5 times the water depth, the double am 
tude is 0.11 times the water depth, and the velocity gradient at mid-dept 
used, the gradient parameter will be 0.10. Use of the velocity gradient at 
depth and all other factors the same gives a gradient parameter of 0.135. 
wave of greater height, that is with double amplitude 0.145 times the water de 
and using the velocity gradient at 0.8 depth the gradient parameter is 0 

If, for simplicity, one assumes that the wind velocity distribution near 
ground may be expressed as the 1/6 power of the elevation, (1) an approxin 
value of the gradient parameter may be computed fora stack. If one consid 
a15-ft diam stack 100 ft high, and bases the gradient parameter on the grad 
at mid-height the gradient parameter will be 0.05. If one considers the lo 
portion where the gradient is steeper, for example, the lower 40 ft, the grad 
parameter would be 0.125. These values indicate that the gradient paramet 
used in the investigation are inthe range to be significant for some enginee: 
problems. 


CONCLUSIONS 


In this exploratory investigation of the velocity gradient effect, test col 
tions were simplified in the following ways: 


1. The ratio of cylinder diameter to water depth was essentially const 

2. The Reynolds number range was chosen so that the conventional ¢ 
coefficient was essentially constant. 

3. The velocity variation was approximately linear giving a constant v: 
for the gradient at all points along the cylinder. 


Based on the results obtained under these restrictions, the following 
clusions appear to be justified: 


1. A longitudinal velocity gradient in the flow approaching a circular cy 
der induces longitudinal currents in the wake zone and along the upstream | 
ment, thus, altering the local drag coefficient at various positions a 
cylinder. 

2. The local drag coefficient is reduced at sections subject to hig 
stream velocity due to the rise in pressure as fluid is supplied into the } 
from the regions of lower stream velocity. 

3. The magnitude of the above effect increases as the magnitude of the 
locity gradient increases. 4 


7 


These exploratory studies have indicated that the velocity gradient ef 
is significant in a number of problems of engineering importance, and 
opened up several questions requiring further investigation. Further s 
anticipated to evaluate the influence of some of the factors that were exc 
in order to simplify this exploratory investigation. 
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STEADY FLOW OF GROUND WATER INTO A SURFACE RESERVOIR® 


By William Haushild! and Gordon Kruse” 


SYNOPSIS 


rediction of the water table position and the amount of water discharged 
é ground water is flowing from an aquifer to a surface reservoir has not 
exact. The nonlinear partial differential equation that describes the shape 
> water table is difficult to solve. 

proximate solutions obtained by two different methods for the nonlinear 
ion are presented. Both approximate solutions agree better with experi- 
al results than does the exact solution of the simplified linear equation 
1€ flow. 


INTRODUCTION 


1e solution of one-dimensional boundary problems in ground water flow 
many applications. For example, the increment of water added to bank 
ge when a surface reservoir is filled, or taken from bank storage when 
ervoir is emptied, can be considerable and the amount is important in a 
mination of the reservoir hydrology. The same solution also applies to 
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certain unsteady cases of ground water flow caused by the digging of dr: 
and to bank storage adjacent to streams that undergo changes in water le 
This paper presents four approximate solutions for the problem. Solutio 
and III are first approximations in which the effect of the unwatering of 
aquifer on the area available for the flow of ground water is neglected. S 
tions II and IV are second approximations that take the effect of unwate: 
into account. 

The solutions presented inthis paper are based upon the following assul 
tions: 

1. The saturated aquifer is of infinite extent and overlays an imperme: 
layer of zero slope. 

2. The Dupuit-Forchheimer assumptions hold. 

3. The change in reservoir water level is instantaneous. 

4. The aquifer is composed of an isotropic and homogeneous material. 


SOLUTION I 


Consider the flow into a reservoir as shown in Fig. 1. The flow throug 
cross section of unit width at the distance x from the origin is 


Q = K (a +i + alse Re Ne 
in which K is the permeability of the aquifer. 
If we let 
i H 
D=dt+ Tio tet oh et See 
then 
zs oh 
Q=K (D + h)= oe Slee ienieme © gutted es me ° 
The equation of continuity for the strip of width dx and time increment d 
3Q aye 
— dx dt-= V <; OU.dxs «00043. eee 
in which V is the specific yield of the aquifer. 
Or, by use of Eq. 3 
32h She on 
KD + 2 h+K(32) = dasa FS a Rania ctrovto soe 
By letting a = x. Eq. 5 becomes 
dx2 D ox t D dx2 ee ee eee 


This differential equation is nonlinear in form. If d is much larger than ] 
can be discarded from Eq. 3. The continuity equation (Eq. 4) then gives 
following linear differential equation: ; 
2 
a 2 
Ox 


_ dh 
ot 


O06 0-18 6 late er” ole ee eet wean 
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Eq. Tis ce the by Carslaw and Jaeger3 as it applies to the flow of heat in 
ids. Ferris* applied Eq. 7 to the case of constantly discharging surface 
ins. Todd solved an equation similar to Eq. 6, by numerical methods. 


Land Surface 


Water Table -at t=0 


>= 


feservoir 


FIG, 1.—GROUND WATER DISCHARGE INTO A BODY OF SURFACE WATER 


} results obtained concurred with data obtained from a Hele-Shaw viscous 
vy model. 
rhe solution of Eq. 7 given by Carslaw and Jaeger is 


H 
hy = - g + Hb eee eee ee cece eees (8) 
hich g, the error function, is 
*/Eat 
a G en Sn eee ere ere 
(i) aa fu (9) 
0 


integral is tabulated elsewhere. 
[his is solution I. When the drawdown is small compared to the thickness 
he aquifer, solution I is an approximate solution of Eq. 6. It satisfies the 


idary and initial conditions: 


x = 0 t>0 hn - -2 
x > © hy = # 
3 er 
; t=, 0 x> 20 hy = 5 


4 
“Conduction of Heat in Solids,” by H. S. Carslaw and J. C. Jaeger, First Edition, 
rd Univ. Press, 1946, pp. 50-51. 

«A& Quantitative Method for Determining Ground-Water Characteristics for Drain- 
Design,” by John G. Ferris, Agricultural Engineering, 31:6, 1950, p. 285, 
“Unsteady Flow in Porous Media by Means of a Hele-Shaw Viscous Fluid Model,” 
avid K. Todd, Transactions, Amer. Geophysical Union, 35:6, December, 1954, 
“Tables of the Error Function and its Derivative,” Natl. Bur. of Standards Applied 


. Series No. 41, October, 1954. 
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SOLUTION I 


A second approximation to the true solution of Eq. 6 may be obtained 
applying the method of Picard.’ This method uses each successive appro? 
mation to approach the correct solution. To obtain solution II, the nonline 
terms of Eq. 6 are computed from the first approximate solution, Eq. 8, a 
substituted into Eq. 6.as known functions. The differential equation, as th 
modified, is then solved, again being subject to the initial and boundary conc 
tions. The process of solution requires that particular integrals be found f 
the known functions. 


q [dhy\2 
A particular integral p;, for the term -=— | ~—] is 


D \ox 
hi 
Py att # 2D Cece. te ee ee 6 
3h 
a ul 
Similarly a particular integral Po» for the term - D hy sal is 
os 
a 1 ohy 
Po BR 2D ae ty Oke ied eee (: 
The expanded form of this particular integral is 
rn x2 a x? 
Bee ee ee ae Oe * eee 
2 2 2D Vr V4at Dt 
x2 
22 ie "a4 at 
i Oo os Ol cee dea es Cees 
Dvr V4at % ( 


An approximate solution of the modified form of differential Eq. 6 that satisf 
the initial and boundary conditions is 


ho = py +: po eg sie ove © Mere (6 teen W tial ® ( 


This is solution II, in which q is the first approximation, Eq. 8, altered to ec 
rect hy for the initial and boundary conditions. 


= H H 2 4H H 
a $b +22) i oe lbee = E a eak 
A plot of Eq. 13 is shown in Fig. 2 with the dimensionless parameter h 
as ordinate and x//4 qt aS abscissa. A complete drawdown is assumed 
that the solution canbe compared with experimental data. The solution for » 
linear case, Eq. 8, is also shown in Fig. 2. 
Keller and Robinson® obtained measurements of the drawdown curve i 
laboratory experiment in a sand filled flume. Water was drained from the si 


7 “Memoire sur la theorie des equations aux derivees partielles et la method 
approximations successives,” by Emile Picard, Liouville’s Jour. de Methematiqu 
Series 4, Vol. 6, 1890, pp. 145-210, . 

8 «Model Study of Interceptor Drains, by Jack Keller and A. R. Robinson, Proce 
ings, ASCE, September, 1959, pp. 25-40. ; ; 


UNSTEADY FLOW 
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All solutions are computed for 
reservoir level lowered to the bot- 
tom: of the aquifer. Ho 


2.0 2.4 


FIG, 2.—-APPROXIMATE SOLUTIONS OF GROUND WATER 
EQUATIONS AND EXPERIMENTAL DATA 
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through a 4-in. perforated pipe located at the floor of the flume. The ent 
flow of water was intercepted by the drain. These conditions simulated f 
into a surface reservoir where the water is drawn down to the bottom of 

aquifer (d = 0). The shape of the water table was determined before the dra 
down curve reached the upstream boundary of the flume. Results were, the: 
fore, the same as for an aquifer of infinite extent. The sand used had vah 
of K = 0.034 fps, and V = 25.5% and was practically homogeneous. The data 
Keller and Robinson are compared with the theoretical approximation, Eq. 

in Fig. 2. 


SOLUTION III 


An alternative second approximation for the bank storage case may be « 
tained by a method of equating flows. A first approximation is obtained as 
solution II and the unit discharge, Q, is computed from it. This flow is tl 
substituted into Eq. 3 and an improved solution is obtained by integrating 
resulting expression. The first approximation will be solution III. It is ¢ 
tained by solving the linearized differential equation 


2th _ dh 7 
yx2 Teche es « ass Cake or lots 
subject to the conditions 
x = 0 pe MH) h = 0 
t = 0 Kae 0 hes 
X + © h-H 


The significance of the notation used in this case is shown in Fig. 3. 
The required first approximation, solution III, is then 


hg syle heen eee Sota een oe A ec ( 


Solutions I and III satisfy the same equations and boundary conditions with o 
the reference point changed. These solutions, therefore, are identical. 


Land Surface 


Water Table at t=O 


Reservoir 


FIG, 3.-GROUND WATER DISCHARGE INTO A BODY OF SURFACE WATER. | 
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SOLUTION IV 


first approximation to the flow Q, through a cross section of unit width 
Ss 


hich D=d+ ts and hg is given by Eq. 15. 
> drawdown is taken into account, the expression for flow is 
e oh 

Q=-Kd@+h > 


btain a closer approximation of h, Eqs. 16 and 17are arbitrarily equated: 


oh ohg 
K (d+h) = KD Spann tee Stee ss (18) 
tegration with respect to x yields 
(4 +h)? 
we ee = D hg + Cy oa ee elie Me tle" ‘ew Were! ts ae (19) 


constant C; may be evaluated for the boundary condition that h = 0 when 
. The value so obtained is 


2 
d 
C; = oie oc ee ee eee ee we ew we we (20) 
, by substitution 
(a+h)? = a? + 2 Bib svete scares ts (21) 


ion IV is obtained from this relation in the form 


hy = y2Dh3 +d? MR EERESS v0 oe ret en ne (22) 


solution satisfies the initial condition h = H when t = 0. If hy is measured 
the same datum as hg so that they may be compared directly, Eq. 21 be- 


Ss 
hy, = Y2Dh, +0" Sat kd 9: Kees ations (23a) 


he case of a complete drawdown it takes the form 


h, =,/Hh, -0.5H = HANG 0.5 He bare: (23b) 
equation is shown on Fig. 2 for comparison with experimental data. 
comparison of these solutions shows that the second approximations II 
V show closer agreement with the experimental data than do the first ap- 
mations I and III. Fortunately, the best agreement is obtained from so- 
. IV, which is the easiest to use. The comparison is made here with the 
me case of complete drawdown. 

l the solutions are approximations based on the eeetnanten that d is very 
compared to H. This assumption is least valid for the case of complete 
lown. Therefore, the results from using this solution with partial draw- 
3 Should be good. 


+ 
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Solutions I through IV were derived for the case of inflow into a reser" 
upon lowering of the reservoir water level. However, the same equations 
solutions apply to the case of outflow if Q is considered to be negative for | 
into the aquifer. 


COMPUTATION OF a 
Curves similar to those of Fig. 2 can be used to compute a for a g: 


aquifer if values of h, H, d, x, andt have been measured in the field. F: 
these values H/D can be computed and curves drawn for solution I or solu 


TABLE 1.—ERROR IN DETERMINATION OF @ 


h/H Solution I, in Solution I 
percent percent 
-0.41 88 84 
-0.12 {i é 36 
+0.06 60 24 
+0.42 21 17 


IV using h/H as ordinate and x/¥4 at as abcissa. For a given h/H, a cal 
computed from the measured value of x/vt and the value of x/V4 at obtai 
from the curve. 

Solutions I and IV were used to determine a bythe above method for seve 
of the points from Keller and Robinson’s data.8 These values were compa 
with the value of a computed from laboratory measurements of K, D, ant 
Table 1 indicates the percentage error in a@ as determined by these two si 
tions. 
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DISCUSSION 
DISCHARGE FORMULA FOR STRAIGHT ALLUVIAL CHANNELS® 


Discussion by R. J. Garde and John L. Bogardi 


R. J. GARDE,! A.M. ASCE.—The authors are to be complimented for their 
rk concerning the discharge formula for the straight alluvial channels. They 
ye shown that if the mean velocity of flow in a straight alluvial channel is 
en by the equation 


n C,, X, and y are not constants but they depend on the size of the bed ma- 
ial and the regime of flow. According to the analysis carried by Messrs. 
land Hwang, C, can vary between 11 and 287, x can vary between 0.35 and 
1, and y can vary between 0.30 and 0.57. For the large volume of flume 
a, the correlation obtained seems to be remarkable. 
During 1956-59, the writer carried out certain investigations? concerning 
velocity distribution and the sediment transport in alluvial channels. One 
he significant findings of these investigations was the fact that velocity dis- 
tution and sediment transport are very much dependent on whether the bed 
the channel is plane, or has ripples and dunes on it. The work of Messrs. 
| and Hwang confirms this finding. 
In order to find out whether this new formula proposed by Messrs. Liu and 
ang can be applied to the naturalchannels and rivers, the writer has carried 
certain preliminary analyses and the results are reported herein. The data 
d for this analysis cover wide ranges of the discharge Q, the mean velocity 
the hydraulic radius R, and the bed material size d. Table A gives a sum- 
ry of these data. 
The procedure followed was: 

Vou Vard 
1. To determine the regime of flow according to ae Se Tes Criterion. 


jote.—This paper is a part of the copyrighted Journal of the Hydraulics Division, 
ceedings of the American Society of Civil Engineers, Vol. 86, No, HY 7, July, 1960. 
‘November, 1959, by H. K. Liu and S. Y, Hwang. 

Lectr. in Civ. Engrg., Univ. of Roorkee, Roorkee (U.P.), India. 

“Total Sediment Transport in Alluvial Channels,” by R.J. Garde, Ph. D. Disserta- 
submitted to Colorado State Univ., Fort Collins, Colo., January, 1959. 
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2. To determine x, y and C,. 
3. To compute the mean velocity of flow by the use of Eq. 1. 


It was found that several points from these field data fell in the region 
‘Transition’ of Fig. 1 (given by the authors). When this predicted regime w 
used and the corresponding values of Cg, X and Y were chosen, the values 
the mean velocity could be computed. It was found: 


a. For points falling in the dune regime, the observed velocity, most of t 
time, was greater than the computed velocity. 

b.. For the points falling intransition regime the observed velocity was le 
than the calculated velocity. 


Considering, for example, the Punjab Canal’s data, 8 points out of 42 fell 
the transition regime and for these 8 points the computed velocities we 


TABLE A.—SUMMARY OF DATA 


Opihi (New Zealand) 
Wanganui (New Zealand)b 
Mountain Creek (U.S.A.)°© 


Canals (U.S,A.)2 


River or Canal Observation] Range of | Range of Range ¢ 
and Country. Points d,in mm, | R, in ft. Q, in ef 
(1) (2) (3) (4) 
_ 

Punjab Canals (India)@ 42 0.15 - Py 
0.43 3.3 

Sind Canals (India)@ 20 0.02- 1.5- 
OVE 3.4 

Ganga Canal (India) 8 0,43 3.4- 
4.0 

Indus River (India) 6 0.97 6.8- 
7.0 

Tiber River (Italy) 11 2.0- 
3,1 
Rakaia (New Zealand) 5.0 
bias 
1.4 
1.5 
2.3 
2.3 
2.6 


ry 


4 «Theory and Design of Stable Channels in Alluvial Materials,” by D. B. Simor 
Ph. D. Dissertation submitted to Colorado State Univ., Fort Collins, Colo., May, 1 F 

b “Slope Discharge Formulae for Alluvial Streams and Rivers,” by E.C. Schnacke 
berg, Proceedings, New Zealand Inst. of Civ. Engrs., Vol. 37, 1951. 

© “Bed Load Transportation in Mountain-Creek,” by H. A. Einstein, U. S. Dept. 
Agric. SCS-TP 55, August, 1944, 


greater than the observed velocities by 35% to 125%. For the 34 points’ 
fell in the dune regime, the computed velocities were always smaller than 


observed velocities by 14% to 43%. Similar tendency was noticed in the c 
of Sind Canal’s data. : 
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Ine of the possible reasons for such high variation (in the case of data fal- 


in ‘transition’ regime) could be that the regime predicted by a - we 


erion is not correct. In fact, the authors have stated that this criterion is 
ed on the flume data for which the Froude numbers are much larger than 
é for natural streams. The accuracy of the results obtained by the Liu- 
ing formula depends on the correct prediction of the regime and therefore 
criterion for determining regime of flow sia di very important factor. 
1S pointed out by the authors, = - _ or a - —- criteria predict the 
mes of flow with fair degree of accuracy for the flume data. However, 
n applied to the natural river data with larger depths and flatter slopes, 
criteria encounter serious difficulties. In the case of natural river data, 


use of high values of shear, W attains a large value and thus the points 


in transition or antidune region. Actually these points belong to ripple 
> region. This discrepency occurs because the Froude number is not taken 
consideration while determining the regime. 

‘rom the point of view of dimensional considerations it can be shown that, 
eet approximation, the ratio of height of dunes to length of dunes depends 


2 ee 


¥gd ve D- 


the Froude number is an important factor in the determination of regimes 


W. B. Langbein® and Z. Matunashi? have also shown 


q is important in the determi- 


Ow in alluvial channels. Similarly, ae 
Ss 
m of velocity distribution in alluvial channels.2 For these reasons one 


d expect that 


a = qa and F,. will be significant parameters for establish- 
a 

the criterion for the determination of regimes of flow. This was done by 

Garde and M. L. Albertson? 5 and this criterion seems to be quite satis- 

ry. It is applicable to flume data as well as natural river data. 

ince the Froude number, which is one of the parameters determining re- 


S of flow, contains the mean velocity of flow, this -.F, criterion 


b 

d 
ot be used while determining the mean velocity of flow. However, this 
Trion can be used with advantage in the present case where the purpose is 
udy the applicability of equation 1 to natural rivers and canals. 


or these reasons, the regimes of flow of all the data were determined ac- 


ing to aT 2 ae F; criterion; several points which fell in the ‘transition’ 
Ys 
Vv . 
ne according to a - we criterion, now fell in the ripple-dune regime. 


made a considerable difference in the computed mean velocity of flow 
rding to Eq. 1. Table B shows the comparison between the observed mean 
ity and the computed mean velocity for the regime predicted by the two 


Hydraulic Criteria for Sand Waves,” by W. B, Langbein, Transactions, Amer. Geo- 
val Union, 1942, pp. 615-618. 

1. Matunashi, Personal letter. 

Characteristics of Bed Forms and Regimes of Flow in Alluvial Channels,” by R. 
‘de and M. L. Albertson, Civ. Engrg. Sect., Report CER 59 RJG 9, Colorado State 
Fort Collins, Colo. n> 
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TABLE B.—COMPARISON BETWEEN OBSERVED AND COMPUTED MEAN VELOC 


0 


—— 


Canal or river Observed V, in fps. Calculated V, in fps 
Transition regime | Ripple-dune reg 
(1) (2) (3) (4) 
Sind Canals 2.41 4.70 2,52 
2.95 5.76 2.83 
= a oo 
Tiber River 3.10 6.50 128 
6.45 92 
6.08 
5.90 
OF 


INNoTATION = 


| PUNJAB CANALS 
INDUS RIVER 
| TIBER RIVER 
RAKAIA RIVER 
OPIHI RIVER 
WANGANU! RIVER 
GANGA CANAL 
MANGLES RIVER 
SIND CANALS 
MOUNTAIN CREEK 
U-S. CANALS 


@#¢04>rRP@00 et 


COMPUTED VELOciITy FPS 


LINE OF PERFECT 
AGREEMENT 


OBSERVED VELOCITy FPs 


FIG, 1.—COMPARISON OF OBSERVED AND COMPUTED MEAN VELOCITIE 
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Thus it can be seen that computed velocities according to ripple-dune a 
in better agreement with observed velocities. Fig. 1 shows the comparisi 
between observed mean velocity of flow and computed mean velocity for tl 
various data. It can be seen that the general tendency is for the observed v 
locity to be greater than the computed mean velocity. From Fig. 1, it can al: 
be seen that 35% of the data have less than 25% error. 

By slight adjustment of the C, - d curve a better agreement could be ol 
tained. The following procedure was adopted. It was assumed that the valu 
of x and y obtained from the authors’ graphs were correct and that adjustme 
could be made for Cg values in such a manner that Observed mean velocity 
flow is equal to the calculated mean velocity of flow. Thus 


Cg computed = BE ag crms Bi ty Sar ( 


R* sv 


New values of Cg were computed in this manner for all the data availab 
to the writer. When there was more than one observation for the same size 
bed material and the same regime of flow, the average new value of Ca wi 
found. Fig. 2 shows the variation of Ca with d for the ripple-dune regim 
Lane and Carson’s data used by Messrs. Liu and Hwang are also plotted 
Fig. 2. Tentatively a mean curve passing through these points is drawn. 1 
to the size of bed material 0.5 mm or 0.6 mm, this new curve is very close 
the one proposed by the authors. However, for coarser material the propos: 
curve gives lower values of Cag than the one given by the authors’ curve. Mo 
data are needed to be plotted on this figure, especially for coarser materi: 
in order to confirm the location of the proposed curve. 


JOHN L. BOGARDI.5—The attempt of the authors to investigate this e& 
tremely involved problem is very much appreciated. The writer agrees wi 
the statement that with the effect of sediment transport upon resistance bei: 
unknown, empirical relationships rather than theoretical solutions should pr 
vide the correct approach. The soundness of this attitude is corroborated | 
the successful study of the authors. Va d 

An almost complete relationship between the three invariants — i 


’ 


K is revealed by Figs. 4 and 7. d 
In these two figures, the points representing identical «< values obvious 


always result from the same experiment, in which the particle size, the sp 
cific gravity of the sediment, and the water temperature are uniform, or 
which there are but slight changes in the water temperature only. Each of t 
parallel straight lines in Figs. 4 and 7, sloping at 2 = 0.555 and Q = 0.56 
respectively, represent the result of an experiment presumably carried c 
under different conditions of observation. Nevertheless, a positive relationsk 
can be established in Figs. 5 and 8 between the coefficient “A” and the invar 


d , 
ant = which corroborates the general validity of the relationship, ter a 


aS empirical correlation by the authors. 


6 Head, Hydr. Lab. of Research Inst. for Water Resources Development, Bud 
Hungary. 


7 DISCUSSION 27 


In the writer’s opinion the explanation given bythe authors asto the compo- 
on of the invariant group K is slightly ambiguous. According to this, the 
uposition of K should be supported by the following considerations: 


l. To conform with the existing knowledge of boundary resistance. 
2. To correlate the data consistently. 
Namely all the variables 


1 the exception of the unknown velocity V are contained in the invariants 
and md 

v 
n principle, there could be no objection to using a variable, respectively 
Tiant, differing essentially from K. Subsequently the authors themselves 
it that the good applicability of the invariant group K (namely that the 


its representing identical ue values should be straight lines) has been at- 


ed by the appropriate selection of values A, m, and N. 
n fact, a wide variety of solutions can be conceived for a function of the 
a 


‘he circumstance, that the above invariants will be involved in the function, 
been illustrated by the authors themselves by the aid of dimensional anal- 
. On the basis of a few measured data the writer has thus temporarily de- 
ined the function of Eq. 1 with fair accuracy by six-variable correlation. 
ypears that with more data available the value of - , respectively of V . 
d be determined graphically with very great accuracy. 

egardless of these considerations, the writer wishes to state that the au- 
s have made very significant advances in this difficult problem. The writer 
es with the authors’ statement concerning the exponents of x and y. In his 
ion the method used for determining x and y is such as to already contain 
ffect of various regimes of movement. It is likely, therefore, that x and 
> related to the invariants that are characteristic for movement. 

he attitude of the authors towards the application of dimensionless analysis 
give rise to misunderstanding. Actually, the use of dimensional analysis 
>does not guarantee that there has not been any important variable omitted 
the investigation. : 


= 
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In connection with the authors’ comments on Fig. 1 in its originaland pre 
ent form, the writer would like to point out, that by representing the relatic 
ship —- 

eA aie a Rates ie ede ee 
Vx 4 


the tractive force can in all probability be determined in an even simpler f or} 


THE SETTLING PROPERTIES OF SUSPENSIONS? 


Discussion by Charles G. Gunnerson 


CHARLES G. GUNNERSON,! F. ASCE.--The author has pointed out that the 
lyses of settling velocities of various sewage effluents shown in his Fig. 2 
‘e made in connection with studies of marine disposal in the vicinity of Los 
reles. During the same period observations were made of the behavior of 
se effluents in the receiving waters. The results of the marine studies have 
npreviously reported.2 It is sufficient here to summarize the data relating 
he sedimentation of effluent suspended solids in sea water. 

The sanitary significance of sedimentation of suspended solids from the ef- 
nt of the Hyperion Treatment Plant lies in the concurrent reduction of en- 
ic bacteria populations in surface waters. Accordingly, observations were 
le of the rate of disappearance of coliform bacteria from the surface. The 
bable effects of dilution and mortality were determined and the effect of 
imentation was computed. 

[it was found that, on the average, coliform bacteria were removed from 
face waters by sedimentation at rates such that 90% removals occurred as 
OWS: 


Hyperion secondary effluent (plus elutriation overflow)—21.0 hr. 

Hyperion primary effluent (plus elutriation overflow)—5.3 hr. 

Orange County Sanitation Districts primary effluent—2.0 to 2.4 hr. (A range 
Se is given because of some uncertainty indetermining the effect of di- 
on). 


se results are based upon samples taken from the upper 1-ft layer of the 
er over periods of up to 1 day. 

n order to compare the field results with the author’s settling tube studies, 
data in his Fig. 2 have been replotted to show suspended solids reductions 
1 function of time by removing the 40 cm depth in the z/t function. It is 
umed that the sedimentation velocities represent removals from a slice at 
40 cm depth. It is further assumed that such removals can be compared, 
east on a gross basis, with removals from the essentially homogeneous 
er 1-ft layer of the ocean. 


December, 1959, by Ronald T. McLaughlin Jr. 

Civ. Engr., Bur. of Sanitation, City of Los Angeles, Calif. 

“Sewage Disposal in Santa Monica Bay, California,” by C. G. Gumnerson, Proceed- 

ASCE, Vol. 84, No. SA 1, February, 1958, pp. 1-28, and Transactions, ASCE, Vol. 
1959, pp. 823-851. 

“Studies on Coliform Bacteria Discharged from theHyperion Outfall, Final Bacter- 
cal Report,” by S. C. Rittenberg, Allan Hancock Foundation, Univ. of Southern 

fornia, Los Angeles, Calif., 1956. 
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Fig. 1 shows both the laboratory and field data for the removal of suspend 
solids with time for Hyperion and Orange County Sanitation Districts effluen 
Although there are minor differences in the groups of effluents tested, 1 
rates of sedimentation in the ocean are clearly greater than those inthe settli 
tubes. Qualitatively, at least, the effect of flocculation of the colloidal fracti 
of the suspended solids appears to account for the greater part of solids r 
moval in the ocean after the first few hours. 
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FIG. 1.-COMPARISON OF SEDIMENTATION RATES IN 
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The flocculation of stream- carried suspended solids discharged into 


water has been studied by marine geologists for many years. K. O. Eme 
considers that the variable composition of source rocks and diagenetic c 
after deposition on the ocean floor are probably more important factors 

distribution of clay minerals in marine sediments off the California coast 


4 «The Sea Off Southern California,” by K. O, Emery, Wiley, 1960. 
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differential flocculation and deposition. Nevertheless, Tj van Andel and 
Postma? have recently used the results of V. G. Whitehouse, et al.,6,7,8,9 
xplaining the near- and off-shore distribution of illite and montmorillonite 
pectively, in the Gulf of Paria. 
Some of the aspects of Whitehouse’s results in studies of differential floc- 
ation and settling of clay minerals appear to be applicable to sedimentation 
ffluent solids in the ocean. He found that the median settling velocities in- 
ased by about two orders of magnitude upon flocculation in saline waters 
n 18 parts per thousand “/gg chlorinity for the various clays tested. Ini- 
Settling rates increase byabout 2% and 22% for kaolinite and illite, respec- 
ly, as the chlorinity is increased from 0.5 to about 3 /00 and are unaf- 
ed by further increases in chlorinity. In contrast, as the chlorinity in- 
ases from 0.5 to 18 “/oo, initial settling rates for montmorillonite are in- 
ased by factors of from 3 to 37 in a temperature range of 26° to 6° C. 
vever, the velocities are still of the order of 1/10 those of the other clays. 
effect of a temperature reduction from 21° to 9° C (this is about the sea- 
al range in Santa Monica Bay in the upper 200 ft of water) is to decrease 
ling rates of the three clays by about 30%. Ions of heavy metals increase 
settling rates; the highest increase of 85% was reported for montmorillo- 
in the presence of 4 x 10-4 molal Fe3+. The presence of carbohydrates 
meses the settling velocity of montmorillonite by 35% and 25% at 0.5 and 
%o0 chlorinity, respectively. Humic acids or proteins cause decreases in 
cities of up to 1% for montmorillonite and up to 30% in kaolinite. As the 
7 concentration decreases, there is a small but significant increase in set- 
¥ velocity. Other factors which influence settling are pH, clay concentra- 
gradients, and the initial state of division of the clay. 
The fate of dispersed organic solids is not as well characterized. Some of 
organics are undoubtedly absorbed onto the inorganic particles, some may 
‘locculated and settled, while the remainder may remain dispersed for a 
siderable time. D. L. Fox, et al.,10,11 have found colloidal or otherwise 
eersed organic detritus in mid- and deep-ocean waters as well as near- 
re waters. It follows that the persistence of this materialin the water col- 
1 is determined by the rate of bacterial decomposition and mineralization 
er than flocculation and settling. 


? 
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The previously cited works are not directly applicable to waste dispos 
Nevertheless, it is reasonable to expect that at least a portion of the efflu 
solids will respond to various environmental factors in manners similar 
those of natural compounds. Certainly, differential flocculation phenome 
must be considered in evaluating the rate and location of deposition of efflu 
suspended solids and their constituents. The converse of this proposition, t 
dispersed, unflocculated solids will be transported for greater distances rm 
have a similar significance. 

One of the conclusions derived from the marine studies of the disappearal 
rates of coliform bacteria from Hyperion discharges was that primary efflu 
was more amenable to final treatment by the ocean than was secondary | 
fluent.2. Two approaches have been used in attempting to verify these resu 
under controlled laboratory conditions. 

In the first, various sea water dilutions of the effluents were placed in la1 
cylindrical columns, and tests were made over a period of time of colifo 
densitites at different depths. It was hoped that the tests would indicate a r 
of sedimentation of the bacteria and a consequent reduction of bacterial po] 
lations at the surface. No such reduction was found, probably because of c« 
vection currents being set up in the column bysmall temperature gradients 
the room and because of the motility of the bacteria themselves. 

The second series of tests involved centrifuging various sea water dilutio 
(The writer is indebted to K. J. Mysels of the University of Southern Califor 
for suggesting this technique). The tests were made at room temperature v 
the centrifuge operating at 1000 g’s for 30 min. Qualitatively, primary efflu 
mixtures showed greater reductions of coliforms than did secondary, thus ci 
firming the receiving water studies. It was also found that greater reducti 
were obtained with higher dilutions; this conforms with the work of Whitehou 
et al., mentioned earlier. 

It was possible to determine the suspended solids reduction with assura1 
only for the 20-to-1 dilution of primary effluent since the absolute solids cx 
centrations in the other dilutions were too small to be accurately measu! 
with standard methods. The solids reduction was about two-thirds that of 
bacterial reduction. 

We now have three converging sets of data, the marine studies, the c 
trifuging experiments, and the settling tube investigations reported by 1 
McLaughlin. The differential settling characteristics of various sewage 
fluents are clearly demonstrated, particularly when effects of flocculatior 
the effluent solids in sea water are included. 

Associated with the utility of centrifuging of sea water, effluent mixtu 
are also indicated. However, it is obvious that this technique requires a g 
deal of experimental and analytical study before quantitative results can be | 
pected that will have the same reliabilityas the very expensive receiving wa 
studies. In any event, it appears desirable to make laboratory tests of sett] 
characteristics of effluent solids which may form bases of design for sew 
treatment plants. Such experiments may well be designed to include the fl 
culating or other effects of the actual receiving water, be they fresh water 
salt water. 
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EARLY HISTORY OF HYDROMETRY IN THE UNITED STATES 
Discussion by Murray Blanchard, Arthur H. Frazier, and J. C. Stevens 


MURRAY BLANCHARD, 73 F. ASCE.—The author has compiled avery inter- 
ng history of hydrometry. Under the heading Great Lakes Investigations 
1aS given a complete resume of the equipment and methods developed and 
i by the engineers of the United States Lake Survey on this work. 

fhe writer had an opportunity to make use of those pioneer developments 
he following hydraulic investigations of the Great Lakes: 


. Discharge from Lake Superior of the St. Marys River at Sault Ste. Marie 
h., for the United States Lake Survey. 

+. Dischargefrom Lake Huron of the St. Clair River at Port Huron, Mich., 
the United States Lake Survey. 

Discharge of the Detroit River between Lakes St. Clair and Erie at 
roit, Mich., for the United States Lake Survey. 

. Discharge of the St. Lawrence River at the footof Lake St. Francis near 
ey Field Ontario, Canada, for the J. G. White Company of New York in 
1ection with Water Power development investigations. 

. Discharge of the Chicago Sanitary District’s Main Channel from Lake 
higan at Lemont, Ill., for The Sanitary District of Chicago in connection 
United States Lake Level Litigation.74 


‘o the hydrometric investigations of the United States Lake Survey referred 
y the author there should be added discharge measurements of the outflow 
ake Superior made in 1901 and 1902 by H. F. Johnson and W. Edward Wilson 
a the International Bridge at Sault Ste. Marie, Mich. 


.RTHUR H. FRAZIER.’5—One finds it hard to realize that this paper on 
srican engineering history has been written by a person who has lived in 
rica for only twelve years. No doubt most of the native-born Americans 
read it found much information therein that was entirely new to them. It 
ccordingly with a humble spirit that this discussion is offered. 

1 only one detail do I find reason to disagree with data contained in the pa- 
The author states that the first discharge measurement made in the State 
laine was made by G. H. Hamlin of Maine State College who measured the 
»bscot River in 1884, whereas George L. Vose of the Massachusetts Insti- 


of Technology wrote’ in 1885 that: 
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“In 1835, several prominent gentlemen in Maine became interested 
in the development of the water power of the Androscoggin River at 
Brunswick; and Mr. Baldwinwas employed to make the necessary meas- 
urements and computations .... From the several gaugings, 4,000 cu- 
bic feet per second were reckoned safe as the discharge which might be 
relied upon as constant through the year for mill power; and with forty 
feet fall in the river at this place, Mr. Baldwin concluded that the Andro- 
scoggin at Brunswick offered a power unsurpassed by any in New England 
at that time occupied.” 


As a matter of fact, it seems quite unlikely that even Baldwin’s measurem 
was the first one made in Maine. A contemporary, and no doubt a close fri 
of his, Ithamar A. Beard, made surveys for the erection of a large cotton f 
tory in Saco, Maine, in 1831. Not long afterward Beard advertised (in 
Lowell Courier during October, 1840) that he, as a civil engineer, offered 
services “in surveying land... ,in measuring and computing the quantity 
water supplied by streams ..., and generally, in any business pertaining 
engineering.” His records of streamflow measurements, and those of simi 
engineers of his time, would provide interesting comparisons with present- 
data if they could be found. It is more than likely that many of Beard’s me 
urements would have been made earlier than any presently on record for 

states in which he lived (Maine, New Hampshire, and Massachusetts). 

In writing a history such as this, one is always faced with the problen 
what items might be omitted for the sake of brevity. That circumstance 
ways provides an almost unlimited area for discussion purposes. Advant 
will be taken thereof in the comments that follow. 

Perhaps “Jean” Nicolas Nicollet was singled out among the explorers of 
Upper Mississippi River because of the magnificent map he made of the a 
at that early data (1836-1439). He was not, however, the first governme 
Sponsored explorer. to have undertaken the job. Colonel Zebulan Pike was. 
dered to search for the head-waters of that river in 1805. General Cass 
later (1832) Henry Rowe Schoolcraft also were sent on that mission. Of th 
three predecessors of Nicollet, only one, Henry Rowe Schoolcraft, actu: 
reached Lake Itaska (the nominal “headwaters” of the Mississippi). It is a 
markable coincidence that a still earlier explorer in this area also bore 
name of Jean Nocolet. In 1634 this earlier Jean Nicolet discovered L 
Michigan, and even undertook to journey up the Fox River in Wisconsil 
quest of the “Father of Waters,” but he did not go far, probably only as far 
Lake Winnebago, and then turned back, only a day’s journey away from his gi 
For the sake of the record, however, a question might be raised regarding 
given name attributed to “Jean” Nicholas Nicollet by the author. Most Ame 
can sources have translated his first name as “Joseph,” rather than “Jez 

From a historical viewpoint it would seem that a report made as of A 
4, 1808 on canals by Albert Gallatin, Secretary of the Treasury, might y 
have been mentioned in the section on such projects, also in the second 
pendix. That report, which reviewed the status of all existing public roads’ 
canals under construction (or completed) as of the date it was written 
available. /7 R 


7 
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17 The Historical Register of the United States, Part I for 1814, edited by T. H.I 
mer, and published by G. Palmer, Philadelphia, 1814. 3 
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nd finally, there seems reason to believe that William Hammond Hall, 
fornia first state engineer, might have deserved a little more attention. 
S mentioned in the section pertaining to the “Progress of Hydrometry in 
West,” but possibly with too little explanation regarding his fine engineer- 
accomplishments. Irrigation received its first recognition on a state-wide 
S in the state of California, and Mr. Hall was the first man to have had 


FIG, 62.—WILLIAM HAMMOND HALL 


ge of that project. So skillfully did he undertake such work that Captain 
Dutton reported. ’8 “In this State (California), irrigationhas made greater 
nces than in any other.” 

S a pioneer in this field, it was necessary for Mr. Hall to develop stream- 
1g methods and equipment as well as irrigation procedures. No doubt he 
the nation’s foremost authority on those subjects during his period of ad- 
stration. An example of his design work on stream-gaging equipment is 
lesign of the first American cable-way for making streamflow measure- 
8. He also improved the design of two current meters, the first of which 
ssistant, C. E. Grunsky described as a “Henry type.” The second current 
r on which he made improvements was of the Small-Haskell type. He 
4d several of this second group to the U.S. Geological Survey in later 
s for conducting its work in Colorado and Nevada. Surely, William Ham- 
1 Hall deserves having his picture in this “Hall of Fame” as much as his 
stant, C. E. Grunsky. 


C. STEVENS, ’9 F. ASCE.—In the writer’s opinion, this paper is the hy- 
lics classic of the 20th century. It is not a treatise on hydrometry and 


$$ ~ 
First Annual Report of the U. S. Irrigation Survey, (later converted into the U. S. 


gical Survey). , 
Cons. Engr. and Partner, Leupold and Stevens Instrument Co., Portland, Ore, 


& 


36 July, 1960 


hydraulics. It is atreatise on menwho made those sciences click in the U 
States since the beginning of its history. And this is not all. Mr. Kolupai 
preparing a world wide bibliography of Hydrometry from the beginning of \ 
ten history, some 6,000 years ago, to the present. Only one with his k 
knowledge of hydraulics and hydrometry and his versatility in languages « 
successfully undertake such a momentous task. 

As an example of the author’s thoroughness, the writer is aware of the 
lowing incident. C. S. Jarvis®9 has stated: “Some of the sub-normal 
flood stages were increased arbitrarily by order of unscrupulous rulers. 
Mr. Kolupaila’s bibliography this nefarious desecration is pinned down > 
“The rulers of Egypt established land taxes called Kharag, correspondi 
the top of the Nile flood in that year. It was in the best interests of the ri 
to have the Wafa (When the Nile reached the 16-cubit (28-foot) stage the sh 
a guardian of the nilometer after the Arabs conquered Egypt, had to proc 
the “Wafa,” a beginning of the annual fast.) proclaimed earlier and the 
levelread higher in order to entitle them to start the tax collection earlie 
to increase the assessment. The scale was doctored for that purpose by c 
of the Caliph Omar: the cubits between 16 and 22 marks were shortens 
one-half. This fraudulent alteration resulted in riots of the populatior 
spoiled the value of the longest record of river observation extending 
1,300 years.” 

The writer is wondering whether the relation of forests to water sup 
could properly be included in the category of hydrometry. At the turn of 
century the forest-water supply relationship was one of the major topic 
the lecture platform and periodicals. The reason for this lies in the fact 
this was the time when the nation was debating, in Congress and the pres: 
proposal that the government should acquire large areas of forest land: 
administer them as National Forest Reserves. The conservationists wer« 
ging their acquisition, but there was one hurdle congress had to make be 
such a course could be legalized. 

Congress had authority over river navigation, but unless it could be s 
that forests were a substantial aid to navigation congress was legally po 
less to acquire and administer forest reserves. The ardent conservatio 
therefore propagandized the nation with all the weapons at the command. 

Two examples will be cited: 


1. Gifford Pinchot, probably the greatest exponent of forestry the vy 
ever knew states:81 “The connection between forests and rivers is like 
between father and son. No forests, no rivers.” From this it would ag 
that forestry is a substantial adjunct to hydrometry. 

2. In the crescendo of direful prophecies, this one from the scrap bo 
September, 1908, wins the blue ribbon: “When our forests are gone, 
streams will dry up, the rivers will cease to run, rain will fall no more 
America will be a desert.” 


One might add “and hydrometry will be but a memory.” One can im: 
Teddy Roosevelt saying to his tennis cabinet: “Go to it, boys, the end jus 
the means. Pull out all the stops. Make America forest conscious. } 
mind the facts. We’ve got to create forest reserves.” 


80 “Flood State Records for the Nile River,” by C. S. Jarvis, Transactions, ; 
Vol. 101, 1936, p. 1026. gS nes 
81 “The Fight for Conservation,” 1909, p. 53. 
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One of the first voices raised against this flood of falsification was a lengthy 
‘scholarly paper82 by General H. M. Chittenden of the U.S. Army Engineers. 
The writer has compiled83 all known long time records of stream flow which 
‘analyzed with reference to the forest cover on their drainage areas. 
Thus, out of a welter of conjecture and falsification, grew one of our most 
uable national assets. National forests now cover one-third of our land 
a. With the sustained yield system now in force not alone by federal but 
0 by many private owners, we are assured a perpetual source of lumber and 
a corollary unrivalled recreation facilities, while soil erosion is kept to a 
limum in those areas. The end did justify the means. 

One tie the forests have with hydrometry is the fact that some forested 
as are actually being cutoff to increase the water supply with the knowledge 
acquiescence of forest management. A forestis a great dissipater of water 
evaporation from limbs and leaves and by the extraction of water from the 
l in growth processes. It is obvious, therefore, if the forest is removed 
re water will be available for storage and use. It appears, therefore, that 
estry in its many phases may definitely be akin to hydrometry. 


ee ae eae 
Transactions, ASCE, Vol. LXII, March, 1909, “Forests and Reservoirs in Relation 
tfansactions: ‘ 2 = 
ream Flow with Particular Reference to Navigable Rivers. 
} Journal of the Assoc. of Engrg. Socs., Boston, July, 1913, “Forests and Their Ef- 
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ERALIZED DISTRIBUTION NETWORK HEAD LOSS CHARACTERISTICS 


Discussion by Paul C. Constant, Jr., Marcel Bitoun, Claud C. Lomax, 
Joseph W. Maier and Thomas C. Miller, and J. M. Robertson 


AUL C. CONSTANT, JR.4—The author describes a mathematical method 
hich pressure conditions of an existing water distribution system can be 
ined easily for different ratios of Qp/ Qq- Before the method can be used, 
> must be available at least one set of measured data for two different 
Yq ratios. Also, the various head losses are calculated under the assump- 
that they are a result of Qp/ Qg ratios where the loads are proportionately 
ged and all other system parameters remain unchanged. 
lthough the restrictions imposed seem quite severe, the author’s method 
alculating different head loss conditions is useful for municipality distri- 
m operations for a “closed” system, that is, the operating end of the busi- 
The author very clearly indicates the usefulness of his method for ob- 
ng intermediate pressure conditions. However, he does not stress that 
cation of his method of calculating head losses is mainly for the operations 
yf distribution problems. An example where the author’s method may be 
effectively is in the determination of pressure gradients in a system for 
ought” year. 
andom checks on the author’s calculations for n and ¢ were made using 
lata from Table 1. Anerror, An, no greater than 5% of n (n = 2.46) was 
|. Further calculations were made to determine the error, Ad, in @ (@ 
344) due to Aninn. This error was no greater than 2%. These calcula- 
indicate that the author’s method is within the limits of engineering ac- 
sy. In fact, the errors are trivial because (a) estimates are made on de- 
is, (b) loads are grouped, and (c) of error in meter readings. 
1€ approximate value of the relative error in (a may be determined 


Q 


y be using the general formula for errors, 
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ich E, represents the relative error; N is the function of several inde- 
mt quantities; jj denotes the parameters (i= 1, 2, 3,..., k); and Ayj 


$e ee 


January, 1960, by M. B. McPherson. j 
Engrg. Dept., Midwest Research Inst., Kansas City, Mi. 


Ph ae 


40 July, 1960 


represents errors in the parameters. The general formula, Eq. 2, can be 


plied to the author’s equation, 
h n 
the a 1) (=) Pa ie fie ah Soe, Sk ye oe a vy ee 
en a 


to give the approximate relative error in the calculation of the function 


1.85 


— \) a+ -[+(2Y oa (#)] a9. 


in which Q, and Qg are treated as constants, and An and Ag are the error 
the calculated values of n and @, respectively. 
st 
h 


The approximate relative error in oes 
d 


was found to be insignifi 
by using Eq. 5, and selecting Ag and An to be the errors previously calculz 

There is an error in the paper under the heading PROPORTIONAL Li 
CHARACTERISTICS WITH EQUALIZING STORAGE where the author so 
for n. The number 0.67 should read 1.67. 

The author could enlarge upon the analysis ofa particular distribution : 
tem with the investigation of flows in the various pipelines. This is a nat 
consequence of having the head losses. The flows can be calculated simpl 
means of the Hazen-Williams hydraulic formula (constant density) in the fc 


nm \1/1.85 
git Be one 4 


in which Q is the flow rate; k, is the head-loss coefficient of the pipelin 
denotes the arbitrary constant a account for flow units; and h represents 
head loss. 

The method of calculations described in the paper are of little value tc 
person concerned with improvements ofa particular distribution system; | 
as, (a) determination of the best locationfor pipelines, (b) comparing alter 
arrangements of systems planned for construction, (c) optimizing a syster 
as to minimize the number of sources required, or (d) selection of pipeline 
ameters for the best combination of economy and performance. These p! 
lems, along with others, seem to be in the majority for consultants in the ; 
of water and gas pipeline distribution analyses. This, in part, was implie 
the author’s statement, “It must be realized that these values of ¢ andnw 
be changed if either the load or the piping in the network of Fig. 5 were 
fied in any way, or if m was changed.” " 


> 
MARCEL BITOUN, ° M. ASCE.—Publication, by the author, of the resu 
the experiment analysis of large, actual networks is very valuable. Such 
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, unfortunately, too scarce. The author is also to be complimented for hav- 
put in evidence definite characteristics of distribution systems with mixed 
es of loading (residential and industrial). 

A lengthy part of the paper has been devoted to the justification of an hy- 
ulic property of proportionally loaded networks. Three distribution sys- 
is of increasing complexity were analyzed and used to empirically justify 
fact that “with proportional loads the percentage distribution of flow in in- 
idual pipes in a balanced network is constant irrespective of the magnitude 
he total demand.” The author expressed this as a “contention” in his con- 
Sions. 

This property has been used by the writer and, it is presumed by other hy- 
ulic engineers for the analysis of balanced networks in which the primary 
rest was in the distribution of discharges rather than the determination of 
d losses. Actual flow values were replaced by proportional figures, the 
zest being assumed to be 100 arbitrary units. Sometimes cumbersome fig- 
Swere thus eliminated. The writer wishes to submit here a demonstration 
he validity of the hydraulic property empirically deduced by Mr. McPherson. 
Consider a distribution system N without equalization storage governed by 
ndary conditions expressed by given ingoing and outgoing discharges. In 
; system, let F be one balanced regime of flow, characterized by a set of 
ding and outgoing discharges, and the corresponding values of flows in the 
ividual pipes that compose the system. 

We have, for any junction of this system, 


i i 
- Q 20 eee (7) 
i=1 


yhich r is the number of pipes connected at the junction. 
n any loop of the system, we have 


n 
Peo dice PO ete errant ake (8) 
i=l 


fhich n is the number of pipes in the loop. 

set us multiply all boundary discharges by a factor q (q # 0), and let F' be 
new regime of flow in the system. 

n any loop Eq. 7 can be multiplied by q and yield 
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Eq. 8 multiplied by q™ yields 


n 

m 
), Kyo hig sso sania: Bee eee ( 
i=1 


If 


then 


Eqs. 10 and 12, valid for any junction and any loop respectively, show tl 


1) the regime F' is hydraulically balanced, 
2) the flow for F' in every individual pipe is equal to the flowfor F in the se 
pipe multiplied by the factor q. 


Moreover, this flow distribution is the only one compatible with the gi 
boundary discharges. Demonstration of the unicity of the flow distributio1 
a system where boundary discharges are given was furnished by Ang 
d’Auriac. 12 

Therefore, since the regime F', in which all discharges are those char 
terizing F multiplied bya constant factor, is a solution, and since this solut 
is the only flow distribution in the network corresponding to the proportic 
load, we have shown that the flow in each individual pipe is a constant fract 
of the total load, which confirms Mr. McPherson’s findings. 

The number of computations made on the analyzer for the constructio1 
Fig. 4 could therefore have been considerably reduced. A single analysis 
each network would have been sufficient to provide one basic flow distributi 

The second comment which is offered here refers to @ andn. For] 
McPherson’s Eq. 2, ‘ 


yh n 
re (2,24) , 


Q 
it was stated that for ~P = 1, ¢ = K. 


However, two field tests were conducted for each location. In fact, only 
test is necessary. The parameter @ is actually the hydraulic resistance of | 
system, and can be computed from the physical characteristics of the pi 
which compose it, provided, of course, that these characteristics (C, ‘5. Dm 
known. A single field test will furnish the second equation necessary for 
determination of n. 

The last comment refers tothe relative influence of non-proportional va 
ations of the industrial loading and of the residential loading in the same n 
work. It is felt by the writer, although this is not a certainty, that deviati 


12 «4 propos de Vunicité de solution dans les problémes de réseaux maillés,” (At 


the unicity of the solution to network problems), La Houille Blanche, No. 3, May-dJt 
1947, Grenoble, France, p. 209. a 
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.a@ proportional loading within a certain range will not result in important 
rences in the distribution of flow in the system. This would be due to a 
pening” of the difference throughout the system. Comments of the author 
d onhis experience onthis point and especially on results of computations 
> on the analyzer might prove conclusive. 


LAUD C. LOMAX,® M. ASCE.—The author is fortunate in having had an 
rtunity to compare actual flow conditions in a pipe network with the results 
ned on a McIlroy Analyzer. The writer has not had that opportunity and 
_present his comment in the light of this lack of experience. 

rom a theoretical standpoint, the technique of proportional load hydraulics 
9e firmly supported. The comparisons that the author makes between the 
roy Analyzer results and the actual flow distributions and the pressures, 
“some discrepancies. The differences may be attributed, in part, to the 
wing: metering limitations, Fluistor selection limitations, and mistakes 
may have occurred in the McIlroy Analyzer solution. The conformity of 
omparisons is effected by the reliability of the field measurements of flow 
ibution, pressure losses, and data such as pipe sizes, lengths, and friction 
rs. The writer is of the opinion that the results presented in this paper 
vell within the range of discrepancy to be expected. 

the light of the writer’s experience, which has been primarily with dis- 
tion networks for small western cities, a word of caution on general ap- 
tion of proportional load hydraulics is in order. In the smaller cities the 
aration of a skeletonized or arterial network is usually not feasible, be- 
e there may be onlyone run of large pipe from a pumping station or res- 
ir, and the remainder of the network is made up of 4 in., 6 in., and 8-in. 
3. The author stated, “This is the smallest distribution district in Phila- 
lia and the only one in which fire flows governed the design of future ar- 
1 piping requirements.” For smaller cities the fire flows are usually 
‘olling. For western areas irrigation loads are significant parts of the 
system load. Consequently, it is the writer’s opinion that application of 
»ortional loads for design conditions of flow networks in small cities is not 
ble. 

1e author has not so stated, but it is implied that the source magnitudes 
be proportional as wellas the load magnitudes. The smaller cities draw 
both reservoir and well sources, and usually all sources are required to 
fy the maximum demand condition, whereas conditions less than maximum 
re pumps for reservoirs that maybe scattered throughout the community. 
would be very difficult to achieve proportionality of source magnitudes in 
tem where reservoir water surface elevations and pump discharges are 
. One would need additional pump capacities at each pumping station, 
1 were proportional to the percentage contributed by that station as well 
oportional to the total load increase. Fire loads, when superimposed on 
aaximum day demand for irrigation, domestic, and industrial usage, do 
\erease the load on the system proportionately since these fire loads are 
d at specific hydrants selected throughout the community. The use of 
tional load hydraulics in the smaller cities may be valuable, if judiciously 
2d, to determine the capacity of the system for flow distributions greater 
the maximum day demands. A clue to the behavior of the system with 
OS eS ee ee 
ssoc,. Hydr. Engr., R. L. Albrook Hydr. Lab., Washington State Univ., Pullman, 
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fire loads can be obtained in this way. However, the pressures and flow « 
tributions inthe immediate vicinity of those hydrants being used may be gre 
in error. 

Having used the McIlroy Analyzer for a considerable number of flow « 
tributions, the writer, and probably others having had close contact with s 
eral McIlroy Analyzer distributions, is impressed by the feeling one gets 
a “healthy” network distribution system. In a well-designed system the on 
sion of a section of pipe, other than the critical section leading to a reser 
or pumping station, usually does not seriously limit the capacity of the syst 
Through errors in wiring or selection of incorrect Fluistors, one accident 
sees this demonstrated on a network, in that it is not really obvious tha 
error has been made until a check procedure discovers this error. On ot 
networks, which are “weak,” the omission of almost any pipe in the netw 
distribution system will reduce its capacity. Proportional load hydraulic 
probably much more applicable to a “healthy” system than a “weak” one. 

With very low flows, conditions may be present wherein laminar flow ex 
in some of the pipes. To apply proportional load hydraulics and increase 
flow to the point where these flows become turbulent, is certainly unwarran 

Referring to the author’s statement on gas flows, it should be noted that 
flow distributions are related to the 1.85 power of the flow when the Fritz 
or Panhandle formulas are used. 

The author is to be complimented for having presented, with supporting | 
dence of its validity, a technique of analysis which, if judiciously applied, 
serve designers of flow distribution systems. It should be obvious that 
does not require a McIlroy Analyzer analysis nor a complete system anal 
to initiate use of this technique. Having pressures at critical points in a « 
tribution network andassuming the system proportionally loaded, one can ' 
evaluate the total capacity of the system. Of course, one can start wii 
McIlroy network analysis as well. 


JOSEPH W. MAIER! and THOMAS C. MILLER,® A.M. ASCE.—The pro} 
tional loading principle, as presented in this paper, does not seem to be 
plicable to gas system networks. The term “proportional load” is descr 
as the typical design assumption that each individual consolidated load va 
or fluctuates about its mean in direct proportion to the total system fluc 
tions. Thus the maximum and minimum hour demands of the average 
maximum day would be taken as fixed multiples, not only of the annual ay 
ages for the district as a whole, but also for each individual local load. It 
been the writers’ experience that a consolidated load fluctuates accordin 
the characteristics of the primary gas load of which it is composed. F 
this load characteristic fluctuation, it can be further concluded that a f 
multiply representing total system demand variation is apparently not ar 
cable to individual consolidated loads in a gas distribution system. 

The composition ofa consolidated gas load will include gas used for on 
more of the following purposes; commercial, industrial, house heating, ci 
ing, water heating, refrigeration, and drying. A consolidated load consis 
of industrial customers will not realize appreciable yearly variation, whi 
consolidated load composed of residential house heating customers will be: 


ee 


7 Asst, Engr. of Distribution Design, Philadelphia Gas Works, Philadelphia, P. 
8 Staff Engr., Philadelphia Gas Works, Philadelphia, Pa. 


DISCUSSION 45 


ted by temperature and wind velocity variations. It can be seen that each 
€ above mentioned gas usages has specific characteristics which when 
dered in a network analysis cannot be represented by a fixed multiply. In 
ion, Seasonal industrial customers are supplied during the non-heating 
yi. 

) illustrate gas load fluctuation, letus use Figs.8 and 9. Fig. 8 is a sche- 
> flow breakdown of a winter day, in which Q,; is the total hourly sendout, 
fers to the industrial customers, Qg indicates the residential non-heating 
mers, and Q3 is the residential heating customers, while Fig. 9, in which 


ot 2M Q: 


|OM 


WINTER 


FIG, 8 


Q, 6M 
SUMMER 
Q: IM 
3M Q; “a 
FIG, 9 


the residential non-heat customers plus seasonal industrial customers, 

refers to the non heat (base load) of residential heating customers, 
ssents the same network on a summer day. Although the summer hour 
mut is 60% of the maximum winter hour sendout, the consolidated loads of 
ummer day (Qj, Qo, Q3)are not 60% fixed multiplies of the corresponding 
— winter loads. 


M. ROBERTSON, 2 F. ASCE.—It seems desirable to point out that the 
ion of pipe network flow and head loss problems with a digital computer 
Prof. of Theoretical and Applied Mechanics, Univ. of Illinois, Urbana, Il. 
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after the Hardy Cross technique has been improperly termed relaxation. 
though many authors use the term relaxation so as to include the iteration 
cess employed in the Cross approach, there is a distinction that shouid be 
served. The error in terminology is not surprising since the introducti 
the iteration technique of obtaining solutions to such problems by H. Cro 
was but a by-product of his work in structural analysis, where he did intro 
the technique (as noted by G. E. Forsythe,11 a relaxation technique was pr 
bly used as early as 1823 by K. F. Gauss for solving the simultaneous e 
tions encountered in surveying) later named relaxation by R. V. South 

The pipe network analysis technique is sometimes termed the Cross-Dc 
method, Since it was at the suggestion of J. J. Doland12 that Cross consid 
the application of the methods of numerical successive approximation to 
work problems. As noted recently by T. G. Chapman, 13 pipe networks mé 
solved either by iteration or by relaxation methods. The relaxation meth 
primarily limited to paper and pencil type calculations, to which it is id 
suited and for which it was originally devised, as so clearly stated by ( 
Forsythe.14 

In spite of the fact that relaxation techniques are applicable to various t 
of fluid flow problems, which we may now attempt to solve with the help of 
speed digital computers, solutions of such problems on digital computers 
rarely, if ever, relaxation solutions. They are merely iteration solutions: 
two terms are not Synonymous. This is because relaxation employs judg 
and reasoning, if not intuition, at every step. These are beyond the usual r 
tality whichis programmed into a digital computer. Sucha machine can me 
plug along, at a very high rate of speed of course, making one simple ca 
lation after another in accord with a rigid prescribed program. Relax: 
techniques, on the other hand, accelerate the iterative convergence of solut 
through various strategems. Thus, computation is first applied to reg 
where the divergence from a proper local result is greatest. In addition, | 
techniques as block and over relaxation are employed where suitable. § 
judgment is required in deciding when and what technique to employ this 
hardly be programmed into a dumb computer. Or, if it can be, the progra: 
just too complicated, too involved. Quite probably nothing is impossible 
computer, but it is better to do a simple computation many times than to t 
complicated one that can achieve the same computation in one cycle. Thus 
Southwell himself has said, it is best to set a computing machine to work 
an iterative process. It is the present writer’s understanding that such is 
deed the case with distribution network analyses. 

The question of the value of the exponent m in the power-law (someti 
incorrectly termed exponential law) head-loss relation is raised by the aut 
Although this exponent may assume values between 1.75 and 2.00, nety 


10 “Analysis of Flow in Networks of Conduits or Conductors,” by H. Cross, Uni 
Illinois, Engrg. Experiment Sta. Bulletin No, 286, November 13, 1936. 
1 “Gauss to Gerling on Relaxation,” by G. E. Forsythe, (translation and comr 
tary on a letter by Gauss dated December 25, 1823) Mathematical Tables and Other 
to Computation, Natl. Research Council, Vol. V, 1951, pp. 255-258, ; a 
12 “Simplified Analysis of Flow in Water Distribution Systems,” by J.J. Doland, 
gineering News-Record, Vol. 117, October 10, 1936, pp. 475-477, 


13 “Relaxation Methods for Pipe Networks,” by T. G. Chapman, Civil Enginee 
and Public Works Review, Vol. 1, November, 1956, pp. 1247-1248, 


14 “What are Relaxation Methods,” by G. E. Forsythe, Chapter 17 of Modern NV 
ematics for Engineers, Ed E. F, Beckenbach, McGraw-Hill Book Co., 1956. ¥ 
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ySis progresses best if a single value is used for any given analysis. 
ainly the author’s proportional load hydraulics method is not applicable 
veral pipes of different m values occur in the system. Network solutions 
Mques may not be limited to fixed-m (versus variable or several m) sys- 
, but the exigencies of the approaches to the solution, such as are being 
issed, often so limit them. Certainly any paper and pencil analysis (plus 
rule or desk calculator and table of powers of numbers) whether merely 
tive or relaxation, can easily employ different m values for each pipe.13 
Mcllroy network analyzer may be used with m other than 1.85. Different 
> could even be employed for different pipes, but this is rather trouble- 
2». Digital computer iteration solutions have been employed for m values 
00 and less, but with the same value for all pipes in any given network. 
Henry and R. B. Peritz!5 note, however, that the use of non-integer val- 
yf the exponent tied up a large partof an IBM 650 machine for the process 
ising the flow rate to this power. A digitalcomputer cannot do this nearly 
imply as one can with an analog computer (such as a slide rule, for in- 
e). This extra effort reduces the number of pipes that can be handed over 
qumber which can be handed when m is 2. 
1€ programming complexities of using different m values for the several 
; of a system are obviously much greater, since a special Q™ subroutine 
i be needed for each value of m to be employed. Either an extensive table 
ch results for a limited set of m values would be needed or some analog 
liary must be developed to which the digital computer may refer quickly 
Q™ is needed. The first method is probably impractical due to the ac- 
time required, while the second is not available. Since the range of flows 
given pipe is usually small, there seems to be no need to consider pipe 
of variable m. It would thus appear that for some time to come our so- 
iS of network systems is limited to a few fixed values of the head loss ex- 
it. Finally, of course, there remains the practical question of whether the 
ent can be determined accurately enough to justify considering anything 
than an average value. 


7 


‘Computer Evaluates Design of Gas Gathering System,” by L. R. Henry and R. B. 
_ Oil and Gas Journal, Vol. 57, No. 17, April, 20, 1959, Also “Evaluation of 
l-Gas System Design with an Electronic Computer,” ASME Paper No. 58-Pet-28 
»ondensed version see Mechanical Engineering, June 1959, pp. 58-59.) 
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THE FOURTH ROOT n-f DIAGRAM@ 
Discussion by Gerald Lacey and N. Rajaratnam 


SRALD LACEY.!1~—The author in writing his paper has, manifestly, two 
ts in mind: first the presentation of a new friction-factor diagram for 
se of practical engineers, and secondly that of challenging the validity of 
ithmic flow formulas. 

-ovided that the engineer makes gooduse of all the information at his dis- 
. it would seem that it makes little difference whether he uses the Moody 
"am which is effectively that of Colebrook and White, or the Blench dia- 
- The author suggests that the latter should be employed and judged by 
ts. There is little doubt that the author’s diagram can be usefully em- 
d but such use can neither prove not disprove the validity of the loga- 
1ic flow formulas which, at present, hold the field. 

10Se who have made a study of recent publications on the theory of flow in 
might well be inclined to conclude that, so far as closed conduits are 
srned, theory is now approaching finality. The use by Colebrook of the 
“Smooth and Rough Pipe Laws” supports this view. It is these laws, so- 
d, in their logarithmic form which the author believes to lack a completely 
al foundation or verification in the range of practical data. Such “proofs,” 
ore properly physical explanations, as have been put forward in the aca- 
¢ world, he regards as patently illogical and mathematically obscured. 
thor puts forward as an alternative his Fourth Root Theory based on the 
late of a universal flow formula for all boundary types which, it would 
ir, would embrace open channels with moving boundaries as well as pipes 
rigid boundaries. The author’s diagram, practical although it may be, is 
tially a vehicle for the presentation of the author’s theories of flow. They 
t lack originality. 

ie author draws attention to the fact that Prandtl was careful to explain 
nis proof of a logarithmic velocity distribution did not apply to pipes but 
ild have been of value if the author had stated when, if ever, the proof 
‘alid—in the case of an open channel for example. The point is relevant 
ich as the author, is support of his fourth root theory as applied to pipes, 
0 hesitation in making use of observations on open irrigation channels, 
porting sediment, and with moving boundaries. 

far as rigid boundaries are concerned it is wellestablished that with full 
lence and low values of the relative roughness the logarithmic equation 
e replaced over a large range by the Manning formula with little error. 


> COO OO TT 


fanuary, 1960, by T. Blench. 
Cons., Sir M. MacDonald and Partners, London, England. 
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The contention of the author that it is “believed” that Manning’s n is sli 
dependent on the size of pipe is debatable. 

When the author refers to “absolute roughness,” and the well known fz 
irrigation canal practice that Manning’s n for exact design has to be rec 
for large channels of the same apparent boundary roughness, he stands on 
ground. The writer more than 30 yr ago (about 1930) was responsible fo 
formula on which the author relies 


1 
ear 


(this equation being in metric units), and commented that 


“the coefficient was one of absolute rugosity, equal to that of Kutter ar 
Manning . . . when the hydraulic mean depth was 1 metre, but remainir 
a true constant so long as... . the size of the silt particle was constan 


The need to vary the rugosity coefficient with the size of the channe 
vanished.” 


This equation however applies to open channels with moving boundaries, 
possibly a very different velocity distribution from that in a channel w: 
rigid boundary. 

There is no doubt that the author’s diagram has the merits of great sim 
city. Whether the practical engineer should be weaned from his predilec 
for the Manning equation when dealing with pipe flow and should use the fo 
root equation instead is an open question. The use of the fourth root throug 
may appeal to a tidy mind. 

When the author states that he is content to treat his value of “e” as a ¢ 
number “recognized intuitively as an observer recognizes faces,” he will 
many who will agree with him. The difficulty in measuring roughness qua 
tatively is a stumbling block and in this matter judgment and experience 
of value. This may be some comfort to the practical engineer but he w 


like some assurance as towhether he should use the author’s “e” or the n 
usual “k,” 


N. RAJARATNAM.13_-Blench is to be highly commended for presenti 
new method of solution for the problems of flow resistance in pipes. 
The writer wants to make a few observations on the projected universa 


of the fourth root n-f diagram. Firstly, the author seems to believe that 
Blasius law, given by the equation 


describes the turbulent flow resistance in smooth pipes even for values 
Reynolds number higher than 10°. But it is widely known that for values 
Reynolds number higher than 105 the flow resistance is better approxin 
by the Prandtl-Karman equation 


12 Proceedings, Inst., Civ. Engrs., Vol. 229, pp. 372, 373. <4 
13 “Research Fellow, Civ. and Hydr., Indian Inst. of Science, Bangalore, India, 
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1 
ays asloge Raff .-0:8 


st, even in the lower range, Eq. 9 agrees with Eq. 8 quite closely!4 so 
when a Single law of resistance is desired, it is Eq. 9 that is normally 
ed. 

condly, regarding the universal adaptability of the fourth root n-f equa- 
the author has suggested that if, in the equation 


oa = 2.0 os Biwi vant de ae 


hickness of the laminar sublayer (which is believed to exist near the 
th boundary for turbulent flow) be substituted for x the Blasius formula 
be obtained. The thickness of the laminar sublayer is givenl5 by the 
ion 


ich 6 is the thickness of the laminar sublayer, u, is the shear velocity, 
is the kinematic viscosity of the fluid. 


ituting 


uy = V 


CO] +) 
2 
5 
5 
— 
ry 
iw) 
— 


ich V is the mean velocity of flow in the Eq. 11. 


d) -Rvt 
ies 11.65 Pei eete PEE aie ee (13) 
ETC i see eas he ae ae (14) 
ituting Eq. 13 in Eq. 11 
R1/4 £5/8 =i CONStANEY & sts te ute is, aedioaitia eis (15) 


Ainst R1/ 4¢= constant which is the Blasius law. 
e derivation!§ of the equation 


= E75 es ache gee sR Oa | . (16) 


e 
lating the values of f for Mannings, Nikuradse and Fourth Root formulas 


= 150 is not clear. Further, the modification of Eq. 16 into the equation 


‘Boundary Layer Theory,” by H. Schlichting, pp. 402, 413. 

‘Elementary Mechanics of Fluids,” by H. Rouse, a) p. 194 b) p. 207. 

Unification of Flow Formulas,” by T. Blench, Proceedings, I.A.H.R., Vol. ie 
957. . Bs 
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seems to rob e of any possible physical significance. Eq. 17 at the « 
looks like one of the many exponential formulae proposed at different t 
For example, in recent times, R. M. Advani!” has proposed a formula th: 


Vf ‘ 


oe © 6/6 © © 65 Bb eee 


in which Kg is (Colebrook’s) equivalent sand roughness. 
Based on Eq. 18, it can be shown that 


8 270 Sue oo oe ae Bale ene 


K, = (31.7 n)® ¢t 


which agrees closely with the formula proposed by the writer for open ch 
flow. 

It is generally known that n decreases with increase in the size of the 
as recently confirmed by the careful experiments of M. J. Webster and 
Metcalf.19 But the variation in n seems tobe very small. It is not clear | 
how the author made n independent of the size of the pipe by the multiplic 
factor d1/12, webster and Metcalf data does not confirm this as shown b 


n (Maximum) d (ft) n ai/12 
0.0246 3.0 0.02695 
0.0244 5.0 0.02790 
0.0236 7.0 0.02770 


The equation for the transition line to rough turbulent flow has been 
as 


It has been shown!5)P that if K is the Nikuradse’s sand roughness and 6 i; 
thickness of the laminar sublayer, the transition starts at = <4.0 and en 


5 > 6.0. Based onthis finding, the equation of the transition line in the M 
Chart has been shown to be 


R vi 
d/2K ~ 400 sitar taelMeNeties aris te ere, 


This equation has been shown to hold good in many cases excepting a few 
crepancies recently reported. Since, in this present case, e€ seems to giy 
17 “Recent Developments in the Flow through Pipes and a New Method of ae 


“a Rough Turbulent Flow through Pipes,” by R. M. Advani, Irrigation and Power, 
er, 1957. 


18 “A Contribution to Turbulent Flow in Open Channels,” by N. S. Govinda Ra 
N. Rajaratnam, Presented at the 30th Research Committee Meeting of the C,B.I, 
India, June, 1960, ; 
19 “Friction Factor in Corrugated Metal Pipe,” by M. J. Webster and LeRe 
calf, Proceedings, ASCE, Vol. 85, No. HY 9, September, 1959, 
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ct magnitude of the roughness projections, the same line of argument may 
9e valid and hence the validity of Eq. 20 is doubtful. 

-is agreed that the Colebrook-White formula, on which the curves in the 
Sition range are based, does not give dependable results for all surfaces. 
ct Colebrook29 himself said that his equationrepresents a mean curve for 
aces of “tar-coated cast iron, wrought iron and galvanised iron.” It has 
claimed by Morris2! that the concept of the “equivalent sand roughness” 
yt sound because Nikuradse’s data on whichit is based gives a rising curve 
€ transition region as against a falling curve predicted by the Colebrook- 
€ equation. Further, the Colebrook-White equation does not make any 
ision for conduit surfaces having either a rising or a horizontal charac- 
stic in the transition region. Thus it is clear that Moody’s Chart does not 
dependable values of f for all surfaces in the transition region. But this 
culty is still there in the fourth root n-f diagram also. Hence a new ra- 
il method of approach should be sought, having an eye for the physical pic- 
of the flow. In this direction, an admirable attempt has recently been 
e by H. M. Morris.21,22 

he writer suggests thatin cases where the flow is of the wake-interference- 
as in the case of many surfaces in open channel flow, it may be advisable 
stain the concept of equivalent sand roughness until a better solution is 
ved. 

egarding problem 6, which was worked out by the author, it is normally 
8 that due to the additional effects of the free surface, the existence of 
ndary currents and non-uniform boundary shear, the equations of resis- 
> for open channel floware likely to be different from those for pipe flow. 
d on the theoretical equations of Garbis H. Keulegan and the experimental 
Its of Mr. Powell, the writer18 has formulated a set of rational formulas 
low resistance in open channels with a new conversion formula for n and 
nd another formula for the transition line to rough turbulent flow. 
cknowledgment.—The writer is grateful to N. S. Govinda Rao for his 
est in the preparation of this discussion. 


Turbulent Flow in Pipes with Particular Reference to the Transition Region Be- 
the Smooth and Rough Pipe Laws,” by C. F. Colebrook, Proceedings, Inst., of 
ngrs., London, 1933, pp. 133-156. ; 
“Design Methods for Flow in Rough Conduits,” by H. M. Morris,” Proceedings, 
, Vol. 85, No. HY 7, July, 1959. ; 

“Flow in Rough Conduits,” by H. M. Morris, Transactions, ASCE, Vol, 120, 1955. 
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INSERVANCY DISTRICTS AS FLOOD CONTROL ORGANIZATIONS@ 


Discussion by Henry J. Tebow 


NRY J. TEBOW, 2 M. ASCE.—The author has presented a refreshing idea 
veloping additional local participation and partnership in solution—of 
ms “in which a Federal interest is predominant,” such as the problems 
or rivers. 

conservancy district approach is an invaluable tool to the consultant 
ymetimes has difficulty, due to special interests and their desires, in 
= acceptance of an over-all plan, which will best fit the district as a 


re should be rewarding work from a conservancy district operation, its 
tants, and the state and federal engineers charged with coordination of 
1g, design, and construction. The conservancy district will make prac- 
the application of a “New Approach to Local Flood Problems, ”% sorely 
to economize on flood control requirements for the nation. 

Colorado Conservancy District Act of May, 1937, specifically paved the 
r the construction of the Colorado-Big Thompson Project, a Federal 
project designed to divert surplus waters of the western side of the 
=ntal Divide to the “water-short” acres of the eastern slope in Colorado. 
1 repayment planning for this Federal project allowed no credit for flood 
|, operation of the reservoirs to maximize the storage of snow-melt 
does provide flood control. The history of events for the project is as 


Initial Settlement and Irrigation in Area 1859 
First State Survey of diversion ($25,000) 1889 
First Bureau of Reclamation investigation 1904 
Reservation for Reclamation use in Rocky 

Mountain National Park Act of Congress 1915 
Consulting Engineer Report 1933 
Public Works Administration Funds 1935 

Conservancy District Organized 1937 
Project authorized (Senate Document 80) 1937 
Initial construction (replacement dam) 1938 
First power produced (To Public Service Company 

of Colorado) 1943 
First supplemental water delivered 1947 
Initial year of major supply to lands 1954 


ee En 
ril, 1960, by Cloyde C. Chambers. 

al. Engr., Design and Constr. Div., Region 7, Bur. of Reclamation, Denver, Colo. 
2w Approach to Local Flood Problems,” by Herbert D. Vogel, Proceedings, 
Jol. 86, No. HY 2, January, 1960, pp. 53-63, 
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Completion for Supplemental Irrigation 1956 
Project Completion (Final Powerplant of six) 1959 
Project Data 
Supplemental irrigation for 720,000 acres. 
Power developed for pumping and sale 183,950 kw. 
Average surplus water diverted 257,000 Acre Feet. 
Average annual power from the water 3/4 billion KWH. 
Annual income to U. S. Treasury approximately $4 million. 
Cost of $159 million—all to be repaid. 
Power costs to be repaid with 3% interest. 
Project water value 1954 $22 of $41 million crop. 


The Colorado Conservancy District Act of 1937 provides for the org 
tion of conservancy districts by any District Court upon petition of a re 
number of property owners. The court appoints a board of directors th 
the power to appoint officers, to acquire and hold property, to appro 
water, to enter into contracts, to levy taxes and assessments, to allot 
and generally to administer the affairs of the District.4 

As a reflection of the indirect benefits accruling tothe entire area in 
in a conservancy district, the act authorizes an assessment on all re 
personal property not to exceed one-half mill during the construction | 
nor to exceed one mill thereafter, except in case of default or deficiency 
an additional one-half mill is permitted.4 | 

The act further provides procedures for allotting water to indivi 
municipalities, and irrigation districts and the collecting of special tax 
the use of such water. The board has the power to contract with the | 
States government for construction, when authorized by an election of pr 
owners. 

The historical record of events for the Colorado-Big Thompson P 
Shows that the Conservancy District Act of 1937 was effective in initiz 
plan of construction just about 50 yr after initial Surveys were started. 

Western pioneers, starting with the Union Colony and early irrigatio. 
Greeley, Colorado (named after Horace Greeley—author of the quote, “Gi 
young man”) overestimated the river water Supply. Late summer and 
year of drought periled the very settlements at the western end of the 
roads. The 50 or more yr without construction on a water plan were cer 
due to lack of a conservancy district law. 

Native American conservatism, and perhaps difficult basic communic 
can be given as the reason that 20 yr lapsed before the Miami Conser 
District idea went west from Ohio. 

Not so with regard to the hydraulics of the MCD dry dams. Ivan E. 
and the fine engineering papers on hydraulics, published by MCD, c 
to the United States Bureau of Reclamation in the 1930’s. The legal va 
a Miami Conservancy District were slower moving westward. 


4 As reported in the Tech. Record of Design and Const., Vol. 1 of 4, Color 
Thompson Project, Bur. of Reclamation, Dept. of the Interior, 1957, > 
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Ni QO79(HY7), 2080(H¥7), 2081(HY7), 2082(HY7), 2083(HY7), 2084(HY7), 2085(HY7), 2086(SA4), 2087 
SA4), 2088(SA4), 2089(SA4), 2090(SA4), 2091(EM3), 2092(EM3), 2093(EM3), 2094(EM3), 2095(EM3),2096 
=M3), 2097(HY7)°, 2098(SA4)°, 2099(EM3)°, 2100(AT3), 2101(AT3), 2102(AT3), 2103(AT3), 2104(AT3), 
405(AT3), 2106(AT3), 2107(AT3), 2108(AT3), 2109(AT3), 2110(AT3), 2111(AT3), 2112(AT3),2113(AT3), 
144(AT3), 2115(AT3), 2116(AT3), 2117(AT3), 2118(AT3), 2119(AT3), 2120(AT3), 2121(AT3), 2122(AT3), 
123(AT3), 2124(AT3), 2125(AT3). 


JUST: 2126(HY8), 2127(HY8), 2128(HY8), 2129(HY8), 2130(PO4), 2131(PO4), 2132(PO4), 2133(PO4), 2134 
sM4), 2135(SM4), 2136(SM4), 2137(SM4), 2138(HY8)C, 2139(PO4)° 2140(SsM4)°. 


TEMBER: 2141(CO2), 2142(CO2),2143(CO2), 2144(HW3), 2145(HW3), 2146(HW3), 2147(HY9), 2148(HY9), 
149(HY9), 2150(HY9), 2151(IR3), 2152(ST7)¢, 2153(IR3), 2154(IR3), 2155(IR3), 2156(IR3), 2157(IR3), 2158 
{RS), 2159(0R3), 2160(IR3), 2151(SA5), 2162(SA5), 2163(ST7), 2164(ST7), 2165(SU1), 2166(SU1), 2167(WWs), 
168(WWw3), 2169(WW3), 2170(WWw3), 2171(WWw3), 2172(WW3), 2173(WW3), 2174(WW3), 2175(Ww3), 2176 
Wws), 2177(WW3), 2178(CO2)°, 2179(IR3)°, 2180(HW3)°, 2181(SA5)°, 2182(HY9)°, 2183(SU1)°, 2184 
Wws)°, 2185(PP2)", 2186(ST7)°, 2187(PP2), 2188(PP2), 


OBER: 2189(AT4), 2190(AT4), 2191(AT4), 2192(AT4), 2193(AT4), 2194(EM4), 2195(EM4), 2196(EM4), 
A97(EM4), 2198(EM4), 2199(EM4), 2200(H¥10), 2201(HYi0), 2202(HY10), 2203(PL3), 2204(PL3), 2205 
PL3), 2206(PO5), 2207(P95), 2208(PO5), 2209(PO5), 2210(SM5), 2211(SM5), 2212(SM5), 2213(SMS5), 2214 
bMS), 2215(SM5), 2216(SM5), 2217(SM5), 2218(ST8), 2219(ST8), 2220(EM4), 2221(ST8), 2222(ST8), 2223 
MTB), 2224(HY10), 2225(HY10), 2226(PO5), 2227(PO5), 2228(PO5), 2229(ST8), 2230(EM4), 2231(EM4), 
P32(AT4)°, 2233(PL3)°, 2234(EM4)C, 2235(HY10)°, 2236(SM5)°, 2237(STS8)°, 2238(PO5)C, 2239(STS), 2240 
L3). 

MBER: 2241(HY¥11), 2242(HY11), 2243(HY11), 2244(HY11), 2245(HY11), 2246(SA6), 2247(SA6), 2248 
+6), Seniea, bie, es 2251(SA6), 2252(SA6), 2253(SA6), 2254(SA6), 2255(SA6), 2256(ST9), 2257(ST9), 
258(ST9), 2259(ST9), 2260(HY11), 2261(ST9)¢, 2262(ST9), 2263(HY¥11), 2264(ST9), 2265(HY11), 2266(SA6), 
#61(SA6), 2268(SA6), 2269(HY11)°, 2270(ST9). 

MBER: 2271(HY12)¢, 2272(CP2), 2273(HW4), 2274(HW4), 2275(HW4), 2276(HW4), 2277(HW4), 2278 
4), 2279(HW4), 2280(HW4), 2281(IR4), 2282(IR4), 2283(IR4), 2284(IR4), 2285(PO6), 2286(PO6), 2287 
POG), 2288(PO6), 2289(PO6), 2290(POS), 2291(PO6), 2292(SM6), 2293(SM6), 2294(SM6), 2295(SM6), 2296 
/Mb), 2297(wwa4), 2298(Ww4), 2299(Ww4), 2300(WW4), 2301(Ww4), 2302(WW4), 2303(WW4), 2304(HW4), 


305(ST10), 2306(CP2), 2307(CP2), 2308(ST10), 2309(CP2), 2310(H¥12), 2311(H¥12), 2312(PO6), 2313(P06), 
s14(ST10), 2315(HY12), 2316(HY12), 2317(HY12), 2318(WW4), 2319(SM6), 2320(SM6), 2321(ST10), 2322 
me , 2323(HWA4)C, 2324(CP2)C, 2325(SM6)C, 2326(WW4)¢, 2327(IR4)C, 2328(PO6), 2329(ST10)°, 2330 
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RY: 2331(EM1), 2332(EM1), 2333(EM1), 2334(EM1), 2335(HY1), 2336(HY1), 2337(EM1), 2338(EM1), 
1), polled Bein 2342(EM1), 2343(SAl), 2344(ST1), 2345(ST1), 2346(ST2), 2347(ST1), 
{1)C, 2349(HY1)¢, 2350(ST1), 2351(ST1), 2352(SA1)°, 2353(ST1)°, 2354(ST1). ee 
RY: 2355(CO1), 2356(CO1), 2357(CO1), 2358(CO1), 2359(CO1), 2360(CO1), 2361(PO1), ; 
2) Steen mete (8 2366(HY2), 2367(SU1), 2368(SM1), 2369(H¥2), 2370(SU1), 2371(HY2), 
2373(SM1), 2374(HY2), 2375(PO1), 2376(HY2), 2377(CO1)®, 2378(SU1), 2379(SU1), 2380(SU1), 
HY2)°, 2382(ST2), 2383(SU1), 2384(ST2), 2385(SU1)°, 2386(SU1), 2387(SU1), 2388(SU1), 2389(SM1), 
ey. 2391(SM1)¢, 2392(P01)°. 


57(PO2)¢, 2458(ST4)C, 2459 
a6 Kin), 2484(OP1), 24 
SA3), 2473(ST5), 2474(SA3) 


), 2547(ST7), 2548(SU2) 
2556(SA4), 2557(SA4), 
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ath 
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